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FOREWORD 


This  study  was  conducted  for  the  U.  S.  Army  Engineer  District, 

St.  Louis  (SLD),  from  January  to  August  1973-  It  was  performed  under  the 
general  supervision  of  Messrs.  W.  G.  Shockley,  Chief,  Mobility  and 
Environmental  Systems  Laboratory,  W.  E.  Grabau,  Chief,  Environmental 
Systems  Division,  and  J.  K.  Stoll,  Chief,  Environmental  Simulation  Branch 
(ESB),  and  MAJ  W.  P.  Emge,  Office  for  Environmental  Studies.  Dr.  V.  E. 
LaGarde,  Project  Manager,  ESB,  was  responsible  for  design  of  the  project, 
development  of  calculational  procedures,  and  development  of  a portion  of 
the  computer  software.  Mr.  S.  J.  Winfrey,  ESB,  was  responsible  for  the 
development  of  the  remaining  computer  software  and  the  operation  of  all 
stages  of  the  calculational  procedures.  This  report  was  prepared  by 
Dr.  LaGarde  and  Mr.  Winfrey. 

Directors  of  WES  during  the  study  and  preparation  of  the  report 
were  BG  E.  D.  Piexotto,  CE,  and  COL  G.  H.  Hilt,  CE.  Technical  Director 
was  Mr.  F.  R.  Brown. 
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CONVERSION  FACTORS,  BRITISH  TO  METRIC  AND  METRIC  TO 
BRITISH  UNITS  OF  MEASUREMENT 


Units  of  measurement  used  in  this  report  can  be  converted  as  follows: 


Multipl- 


To  Obtain 


British  to  Metric 

feet 

0 . 3048 

meters 

miles  (U.  S.  statute) 

1.6093 

kilometers 

acres 

0.4047 

hectares 

acre-feet 

1232.75 

cubic  meters 

Metric  to  British 

meters 

3.2808 

feet 

SUMMARY 


Several  geometric  characteristics  of  water-basin  regimes,  includ- 
ing basin  size  and  shape,  area  underwater  at  specific  water  depth,  and 
cross-sectional  area,  are  commonly  associated  with  benthic,  plankton, 
and  fish  community  population  structures,  although  little  quantitative 
data  are  available  to  support  the  association.  This  two-volume  report 
describes  a general  procedure  that  was  developed  to  calculate  values  of 
selected  parameters  used  to  define  the  above-mentioned  geometric  char- 
acteristics of  any  water-basin  regime.  The  procedure  was  successfully 
applied  to  yield  quantitative  information  for  those  parameters  for  l8 
side  channels  of  the  Middle  Mississippi  River.  Which  of  the  parameters 
selected  as  quantitative  descriptors  of  the  characteristics  are  best 
indicators  of  animal  community  population  structures  is  expected  to  be 
determined  as  a result  of  other  projects  currently  under  way  at  the 
U.  S.  Army  Engineer  Waterways  Experiment  Station. 

Volume  I contains  a description  of  the  procedure  and  the  results 
of  implementing  it.  Volume  II  contains  a set  of  computer -plotted  con- 
tour maps  for  the  18  side  channels. 
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COMPUTER-CALCULATED  GEOMETRIC  CHARACTERISTICS  OF 
MIDDLE  MISSISSIPPI  RIVER  SIDE  CHANNELS 


VOLUME  I:  PROCEDURES  AND  RESULTS 


PART  I:  INTRODUCTION 

Background 

1.  The  continued  implementation  of  engineering  structures,  such 
as  wing  dams,  is  under  study  as  an  aid  in  maintaining  the  9-ft*-deep 
by  300-ft-wide  channel  between  St.  Louis,  Missouri,  and  Cairo,  Illinois. 
Dams  constrict  the  river  flow  to  midchannel,  so  that  the  scouring  action 
of  the  swifter  midstream  flow  helps  maintain  the  desired  channel  depth. 
The  materials  scoured  from  the  main  river  channel  and  other  solids 
already  in  suspension  are  transported  downstream,  and  portions  are 
deposited  in  slack-water  areas.  The  deposit  of  those  materials  along 
the  riverbanks  has  an  effect  upon  the  aquatic  regimes  along  the  river, 
particularly  side  channels  that  are  open  to  the  river  and  have  flow 
through  them  at  normal  river  stages.  Since  the  side  channels  are  impor- 
tant as  fish  and  wildlife  habitats,  particularly  as  sport  and  commercial 
fish  spawning  areas,  there  is  a need  to  examine  the  scouring  action  and 
transport  within  the  main  channel  and  the  probable  effects  on  geometric, 
chemical,  and  biological  characteristics  of  the  side  channels.  An 
analysis  of  the  relation  between  those  characteristics  and  animal  popu- 
lations can  be  used  as  a basis  for  predicting  possible  changes  in  animal 
populations  caused  by  changes  in  side-channel  characteristics  and, 
indirectly,  the  changes  in  animal  populations  caused  by  maintaining  the 
9-ft-deep  channel. 


* A table  of  factors  for  converting  British  units  of  measurement  to 
metric  units,  and  metric  units  to  British  units,  is  presented  on 
page  ix. 
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Purpose 

2.  The  purpose  of  the  research  program  under  which  this  study  was 
performed  was  to  provide  reference  material  to  the  U.  S.  Army  Engineer 
District,  St.  Louis  (SLD),  for  preparation  of  an  Environmental  Impact 
Statement  relative  to  the  development  and  maintenance  of  a 9-ft-deep 
navigation  channel  in  the  Middle  Mississippi  River.*  The  primary 
purpose  of  the  study  reported  herein,  which  is  one  of  several  conducted 
by  the  U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES)  in  support 
of  the  program,  was  to  provide  a portion  of  the  environmental  inventory 
information  needed  to  analyze  relations  between  animal  populations  and 
the  geometric,  chemical,  and  biological  elements  of  the  project  side- 
channel  areas.  A secondary  purpose  was  to  establish  a comprehensive 
computer-accessible  data  base  containing  Middle  Mississippi  River  side- 
channel  data  for  use  by  analysts  to  calculate  parameter  values  for 
additional  side-channel  geometric  characteristics. 

Scope 

3.  This  study  encompassed  the  definition,  development,  and  use  of 
a general  procedure  to  quantify  geometric  characteristics  of  18  Middle 
Mississippi  River  side  channels  from  St.  Louis,  Missouri,  to  Cairo, 
Illinois. 

Approach 

U.  Several  geometric  characteristics  of  water-basin  regimes, 
including  basin  size  and  shape,  area  underwater  at  specific  water  depth, 
and  cross-sectional  area,  have  been  associated  with  benthic,  plankton. 


* The  results  of  the  program  will  be  published  by  the  Waterways  Experi- 
ment Station  in  a report  entitled  "Environmental  Analysis  and  Assess- 
ment of  the  Mississippi  River  9-ft  Channel  Project  Between  St.  Louis, 
Mo.,  and  Cairo,  111.,"  by  W.  P.  Emge,  and  others. 


and  fish  community  population  structures, * typically  through  the  food 
chain.  For  example,  penetration  of  sunlight  to  the  tasin  bottom  has  a 
major  influence  on  the  type  and  mass  of  aquatic  vegetation,  which,  in 
turn,  provides  hiding  places  for  some  fish  species  and  a major  link 
in  the  food  chain  of  all  fishes. 

5.  The  following  parameters  were  selected  to  describe  quanti- 
tatively the  geometric  characteristics  of  side  channels, 
a.  Center-line  length. 
b_.  Average  width  between  high  banks . 
c_.  Water  volume  as  a function  of  water  elevation. 
d_.  Shoreline  length  as  a function  of  water  elevation. 
e_.  Water  surface  area  as  a function  of  water  elevation. 
f_.  Shoreline  development  as  a function  of  water  elevation. 

g.  Rate  of  change  of  water  surface  area  with  respect  to  water 
elevation  (derivative  of  water  surface  area  with  respect 
to  water  elevation)  as  a function  of  water  elevation. 

h_.  Ratio  of  water  surface  area  to  volume  as  a function  of 
water  elevation. 

i^.  Ratio  of  shoreline  length  to  water  surface  area  as  a 
function  of  water  elevation. 

j_.  Bottom  surface  area  underwater  as  a function  of  water 
elevation  and  water  depth. 

k.  Water  cross-sectional  area  as  a function  of  water  elevation 
at  selected  sampling  locations  (stations). 

Since  values  of  side-channel  characteristics  change  as  a function  of 
water  elevation  within  the  side  channel  (and  therefore  as  a function  of 
season),  the  appropriate  parameters  were  calculated  as  a function  of 
water  elevation  as  indicated  above.  The  calculation  of  parameters  as  a 
function  of  water  elevation  makes  possible  the  analysis  of  the  data 
during  the  major  phases  of  animal  life  cycles;  for  example,  fish 

* C.  E.  Warren  and  P.  Doudoroff,  Biology  and  Water  Pollution  Control , 
Saunders,  Philadelphia,  Pa.,  1971- 

G.  E.  Hall,  ed.,  "Reservoir  Fisheries  and  Limnology,"  Special  Publi- 
cation No.  8,  1971,  American  Fisheries  Society,  Washington,  D.  C. 

C.  D.  Sculthorpe,  The  Biology  of  Aquatic  Vascular  Plants,  Edward 
Arnold,  London,  1967- 
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spawning,  juvenile,  and  mature  phases,  which  occur  during  different 
seasons  of  the  year  and,  thus,  usually  at  different  water  levels. 

6.  The  Middle  Mississippi  River  side  channels  identified  during 
the  research  program  are  indicated  in  fig.  1 and  listed  in  table  1 in 
decreasing  order  of  river  mile  location.  (Side  channels  are  identified 
and  numbered  identically  with  those  described  in  the  report  to  be 
published  containing  the  research  program  results.)  River  mile  locations 
of  the  extremities  of  each  side  channel  and  sampling  locations  (stations) 
within  each  also  are  noted  in  table  1.  The  parameters  listed  in  para- 
graph 5 were  calculated  in  this  study  for  all  side  channels  except  side 
channels  5,  7,  13,  14 , 15,  24,  and  25. 
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Fig.  1.  Index  map  of  Middle  Mississippi  River  side  channels 


7.  It  was  necessary  to  know  the  topography  of  a side  channel 
before  parameter  values  could  be  calculated  for  it.  Since  a direct 
investigation  of  the  topography  of  a water-filled  side  channel  was 
impractical,  a procedure  was  developed  to  simulate  the  side-channel 
topography  mathematically  using  assorted  available  data.  However,  the 
data  had  to  be  processed  prior  to  use  in  the  topography  simulation.  The 
simulation  procedure  described  below  was  used  for  the  18  side  channels 
studied. 

Processing  of  Available  Data 

8.  The  available  data  for  the  side  channels  consisted  of  complete 
aerial  photography  coverage,  limited  fathometer  data,  and  sparse  field 
survey  data. 

Preparation  of  data 

9.  Aerial  photographs.  Aerial  photography  coverage,  made  avail- 
able by  the  SLD  and  used  in  this  study,  consisted  of  the  following: 


Type 

Date 

Scale 

Coverage 

Black- and- wh it  e 
infrared 

Jan  70 

1:12,000 

Main  river  channel  and 
floodplain  in  immediate 
vicinity  of  main  channel 

False-color 

infrared 

Aug  71 

1:12,000 

Main  river  channel  and 
floodplain  in  immediate 
vicinity  of  main  channel 

Panchromatic 

(black  and  white) 

Jun  69 

1:2U,000 

Floodplain 

The  original  intent 

of  the  SLD 

in  having  the 

black-and-white  infrared 

photo  coverage  flown  did  not  include  specific  coverage  of  river  side 
channels.  Since  most  side  channels  are  adjacent  to  the  main  channel, 
however,  the  majority  of  side  channels  included  in  this  study  were 
contained  in  totality  in  those  photos.  Where  they  were  not,  the 
panchromatic  photography  was  used  to  supplement  the  black-and-white 
infrared  coverage.  The  supplementary  panchromatic  coverage  was  enlarged 
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by  a camera  process  to  the  scale  of  the  black-and-white  infrared 
coverage. 

10.  Only  one  of  the  18  side  channels  was  totally  contained  on  a 
single  black-and-white  infrared  aerial  photo.  A complete  picture  of  the 
remaining  side  channels  was  produced  by  constructing  photomosaics  with 
the  black-and-white  infrared  photos,  supplemented  by  the  panchromatic 
photos  where  necessary.  According  to  SLD,  the  photos  were  rectified  and 
were  at  a scale  of  1:12,000,  and  these  statements  were  taken  as  the 
basis  for  assuming  nondistortion  of  the  side  channels'  spatial  config- 
urations on  the  photos  and  for  calculating  the  side  channels'  horizontal 
dimensions  from  the  photos,  respectively.  No  control  points  were  avail- 
able on  the  photos  to  validate  these  statements.  Fig.  2 is  a repro- 
duction of  that  portion  of  the  photomosaic  containing  side  channel 
Osborne  and  a portion  of  the  main  channel.  The  photomosaic  was  made  up 
of  two  black-and-white  inf rared ■ photos , and  is  approximately  nine-tenths 
the  scale  of  the  original.  In  fig.  2 the  notations  in  the  area  of  the 
side  channel  include  locations  where  fathometer  and  survey  data  were 
taken  (paragraphs  16-19)  and  water  levels  at  the  times  of  the  aerial 
photo  coverage  and  the  fathometer  runs. 

11.  The  black-and-white  infrared  coverage  was  used  as  the  basis 
for  side-channel  horizontal  dimensions  data  for  several  reasons.  It 
contained  the  most  complete  coverage,  with  a late  flyover  date,  at  the 
largest  scale  available.  In  addition,  the  coverage  was  flown  during  the 
time  when  deciduous  tree  foliage  was  at  a minimum  and  the  water  stage 
was  low.  Also,  the  black-and-white  infrared  film  energy  response  func- 
tions are  such  that  the  water-land  interface  is  most  easily  interpreted 
on  that  type  of  film. 

12.  The  false-color  photography  was  examined  in  a limited  number 
of  situations  to  resolve  ambiguities  on  the  black-and-white  photography. 
It  was  not  found  useful  for  supplementing  the  data  to  any  extent  because 
the  foliage  was  in  full  bloom  in  that  photography,  making  interpretation 
of  high  bank  difficult,  and  the  water  level  was  higher  than  in  the 
black-and-white  coverage.  In  addition,  false-color  photography  does 
not  always  show  the  land-water  interface  distinctly;  it  is  frequently 
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Aerial  photomosaic  of  side  channel  Osborne 


difficult  to  determine  whether  a flat  sand  structure  is  above  or  imme- 
diately below  the  water  surface. 


13-  The  river  stage  in  the  immediate  locale  of  each  side  channel 
was  calculated  jointly  by  SLD  and  WES  for  the  time  of  photo  coverage. 

Most  side  channels  contained  water  at  the  same  level  as  the  main  river 
channel,  providing  a ready  source  of  elevation  data  within  the  side 
channels  at  the  land-water  interface  locations.  Procedures  described  in 
paragraph  28  were  used  to  deduce  the  elevation  of  pools  not  connected  to 
the  main  river  channel. 

ll*.  The  spatial  extent  of  each  side  channel  was  delineated  on 
transparent  overlays  from  the  aerial  photos  by  photo  interpreters . The 
high-bank  position,  as  indicated  by  mature  willow  growth,  was  used  to 
define  the  side-channel  edges.  The  water  entrance  and  egress  locations 
(ends  of  the  side  channel)  were  defined  at  their  junctures  with  the  main 
river  channel,  or  preferably  at  the  location  of  water- flow  control 
structures  within  the  side  channel  close  to  that  juncture.  An  example 
of  such  a transparent  overlay  is  shown  in  fig.  3,  which  is  a reproduction 
of  the  overlay  for  side-channel  Osborne  at  the  same  scale  as  fig.  2.  It 
shows  the  total  outline  of  that  side  channel,  the  locations  where 
fathometer  data  were  obtained,  the  locations  of  dikes,  the  waterline 
within  the  side  channel,  and  stations  where  the  cross-sectional  area  of 
the  side  channel  was  to  be  calculated. 

15-  Since  the  data  on  the  transparent  overlays  were  at  an  incon- 
veniently small  scale  for  subsequent  operations , the  overlays  were  digi- 
tized and  input  to  the  WES  computer,  and  expanded  overlays  were  computer- 
plotted  at  a much  larger  scale  to  yield  working  base  maps.  Digitizing 
was  performed  through  the  use  of  a manually  operated  device,  which 
recorded  the  two-dimens ional  Cartesian  coordinates  of  points  and  lines 
at  positions  dictated  by  the  equipment  user.  Information  other  than 
Cartesian  coordinates,  e.g.  the  elevation  of  a specific  point,  was 
recorded  through  a keyboard.  All  recorded  data  were  automatically  placed 
on  magnetic  tape  during  the  digitizing  process,  and  were  immediately 
available  for  input  to  computer  software  at  completion  of  the  digitizing 
process.  A computer  graphics  rather  than  camera  expansion  procedure  was 
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used  because  of  the  overall  time  and  funds  savings  and  the  fact  that 
perfect  linearity  and  scale  control  were  assured  in  the  computerized 


expansion  process. 

16.  Fathometer  and  survey  data.  During  May  1972,  personnel  of 
the  SLD  operated  the  SLD's  fathometer  equipment  to  obtain  bottom 
profiles  within  the  side  channels.  Bottom  profiles  for  any  side  channel 
consisted  of  both  selected  profiles  across  (cross  profiles  numbered  1-12 
in  fig.  3)  and  a single  profile  down  the  estimated  position  of  the 
thalweg.  Attempts  were  made  to  locate  the  cross  profiles  at  the  side 
channel  ends  (often  at  or  near  control  structures),  on  both  sides  of 
any  known  underwater  control  structures,  and  at  one  position  between 
control  structures.  Since  the  operation  occurred  during  high  water,  the 
fathometer  was  carried  over  the  maximum  possible  extent  of  the  side- 
channel  width  on  cross  profiles  and  over  all  dikes  along  the  thalweg 
profile.  Because  of  the  high  water,  the  boat  carrying  the  fathometer 
was  able  to  reach  the  high-bank  positions  and,  in  a few  cases,  to 
actually  pass  over  the  high  banks  on  cross  profiles.  The  distances  of 
the  starting  and  ending  points  for  a cross  profile  from  local  reference 
points  (e.g.  willow  line)  were  noted  on  the  fathometer  strip  charts. 

The  inability  to  measure  cross  profiles  over  the  entire  distance  from 
high  bank  to  high  bank  for  most  of  the  profiles  was  primarily  due  to 
trees  protruding  from  the  water  over  the  high  banks  and,  in  a few  cases, 
to  bluffs  forming  the  high  banks. 

17.  A survey  team,  operating  in  August  1972  during  a period  of 
relatively  low  water,  supplemented  the  data.  The  team  located  the 
fathometer-profile  positions  and  took  profile  data  over  the  ends  of  the 
profiles  up  to  the  high-bank  position,  and  measured  the  horizontal 
distance  from  high  bank  to  high  bank  along  the  profiles , the  distances 
along  the  side  channel  between  cross  profiles,  and  dike  elevations.  All 
measurements  were  performed  with  rod  and  transit . 

18.  Figs.  U and  5 are  examples  of  the  profile  data  available  for 
this  study  after  the  fathometer  and  survey  data  were  spliced.  The  fig- 
ures consist  of  one  of  the  cross  profiles  and  the  thalweg  profile, 
respectively,  for  side  channel  Osborne.  The  location  of  the  cross 
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profile  in  the  side  channel  corresponds  to  the  location  of  the  line  in 
fig.  2 where  the  notation  3-3  appears,  and  in  fig.  3 where  the  line 
labeled  3 appears.  The  location  of  the  thalweg  is  not  shown  in 
figs.  2 and  3;  the  procedure  for  locating  it  is  described  in  para- 
graph 2h . The  orientation  of  fig.  5 is  opposite  to  that  of  figs.  2 
and  3 because  the  fathometer-survey  team  traveled  downstream  to  obtain 
the  thalweg  profile. 

19.  Submerged  vegetation  and  debris,  bubbles,  and  equipment  noise 
frequently  cause  high-frequency  signals  to  be  imposed  on  the  true  bottom 
signal.  The  portions  of  the  profiles  derived  from  the  fathometer  data 
(e.g.  solid  line  in  fig.  b and  all  of  fig.  5)  were  smoothed  to  remove 
these  effects. 


Integration  of  data 

20.  The  data  prepared  as  described  in  paragraphs  9-19?  except 
fathometer-profile  elevations,  were  integrated  for  each  side  channel  in 
an  identical  fashion,  resulting  in  a working  base  map  of  the  side 
channel.  This  map,  together  with  the  fathometer-survey  profile  eleva- 
tions was  used  to  calculate  the  topography  data  for  the  side  channel. 

21.  The  base  map  was  formed  by  superimposing  all  available  data, 
except  the  fathometer-survey  profile  elevations,  on  the  computer-plotted 
expanded  overlay  of  the  side  channel  (see  paragraph  15 ) • The  locations 
and  elevations  of  all  dikes  were  noted,  as  well  as  the  locations  of  all 
fathometer-survey  profiles  and  land-water  interfaces , particularly  where 
the  water  within  the  side  channel  was  connected  to  the  main  channel. 

22.  The  positions  of  the  dikes  were  checked  by  comparison  with 
the  dike  structures  visible  in  the  aerial  photography.  The  survey  eleva- 
tion data  for  the  dikes  were  checked  by  comparing  water  elevations  in 
the  side  channels  with  the  dike  elevations  to  determine  whether  dikes 
that  should  have  been  above  water  level  actually  were. 

23.  The  positions  of  the  cross  profiles  were  checked  by  measuring 
the  distances  from  high  banks  to  high  banks  from  the  expanded  computer 
plot  and  comparing  those  distances  with  the  distances  measured  by  the 
survey  team  for  the  cross  profiles  (see  fig.  M.  In  addition,  the 
height  of  one  portion  of  the  profile  relative  to  another  (the  profile 
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shape)  could  be  deduced  by  examining  the  positions  of  land  and  water 
along  the  profile  on  the  aerial  photo.  Deduced  features  included  the 
general  position  of  the  thalweg  relative  to  the  center  of  the  side 
channel  and  the  locations  where  the  profile  should  be  high  (e.g.  above 
water)  or  low  (e.g.  below  water).  The  deduced  shape  of  a profile  from 
the  base  map  was  checked  against  the  fathometer-survey  profiles  as  a 
second  check  to  ensure  that  the  fathometer-survey  profile  was  properly 
located  and  oriented.  Various  data  errors  were  encountered  in  this  data 
editing  procedure.  In  several  instances,  the  orientations  of  profiles 
were  opposite  to  what  they  should  have  been;  the  confusion  in  alignment 
was  easily  detected  and  corrected.  In  a few  instances,  the  survey- 
measured  bank-to-bank  distance  on  a fathometer-survey  profile  did  not 
agree  with  the  scaled-off  aerial  photo  distance.  Since  the  boat  crew 
performing  the  fathometer  measurements  attempted  to  keep  the  boat  speed 
constant  during  a profile  measurement  and  at  the  same  speed  from  one 
profile  to  the  next , the  length  of  the  profile  in  question  was  estimated 
roughly  from  the  horizontal  length  of  the  fathometer  chart  of  the  profile 
relative  to  that  of  neighboring  profiles  for  which  the  horizontal  ground 
distance  was  known.  In  a few  cases  where  errors  were  detected,  the 
survey -measured  distance  was  found  to  be  in  obvious  error  (e.g.  factor 
of  two  incorrect).  In  all  situations  where  data  corrections  were 
performed,  the  photo-interpreted  distance  was  given  preference,  since 
the  chance  of  gross  error  was  less  than  for  the  field  measurement. 

2b.  The  thalweg  profile  data  (see  fig.  5)  contained  notation  as 
to  where  it  passed  over  cross  profiles.  Locations  in  the  profile  where 
the  fathometer  had  passed  over  dikes  were  obvious.  The  dike  elevation 
survey  data  were  checked  against  the  profile  depth  data  converted  to 
elevation  as  one  consistency  check,  and  no  problems  were  encountered. 

The  water  depth  on  the  thalweg  profile  was  checked  against  the  water 
depth  on  the  cross  profile  for  each  cross  profile  where  they  intersected. 
The  position  of  the  intersection  of  the  two  was  chosen  so  that  the 
thalweg  profile  was  toward  the  center  of  the  side  channel,  since  the 
boat  crew  had  attempted  to  keep  to  the  center , and  the  water  depth  at  the 
point  of  intersection  of  cross  and  thalweg  profiles  was  the  same.  Where 
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possible,  the  structure  of  the  side-channel  bottom,  as  seen  on  the  aerial 
photo,  was  correlated  with  the  shape  of  the  thalweg  profile  as  a counter- 
check of  the  data  (in  particular,  a check  on  the  correct  horizontal 
positioning  of  the  thalweg  profile).  No  extensive  checking  of  this  type 
could  be  performed,  however,  as  the  thalweg  profile  was  typically 
measured  along  the  deepest  region  of  the  side  channel,  and,  therefore, 
the  majority  of  the  profile  position  was  underwater  even  in  the  aerial 
photos  taken  during  low  water.  No  particular  problems  were  encountered 
where  correlations  could  be  made. 

25.  There  was  a time  lapse  of  28  months  between  the  aerial 
photography  and  the  fathometer  measurements . As  a result  of  the  data 
editing  process,  it  was  discovered  that,  insofar  as  it  could  be  deter- 
mined, the  two  data  sources  were  in  close  agreement  as  to  the  side- 
channel  bottom  structure.  Any  differences  in  shoaling  patterns  between 
the  two  measurements  were  not  extensive.  Channels  did  not  move  from  one 
side  to  the  other  in  the  side  channels  , nor  did  the  shoaling  patterns 
appear  to  change  extensively. 

26 . In  addition  to  including  on  the  base  map  the  locations  and 
elevations  of  all  land-water  interfaces  where  the  water  within  the  side 
channel  was  connected  to  the  main  channel,  so  that  the  water  elevation 
was  known,  various  elevation  data  were  interpolated  onto  the  map  using 
a combination  of  the  measured  and  interpreted  data  as  a basis.  This 
step  was  necessary  to  make  the  full  informational  content  of  the  measured 
and  interpreted  data  available  to  the  computerized  topography  calcula- 
tion. The  interpolated  data  consisted  of  high-bank  elevations  and 
elevations  of  the  land-water  interface  for  pools  within  the  side  channel 
not  connected  to  the  main  channel. 

27.  Measured  data  on  the  high  banks  consisted  of  field  survey 
measurements  on  the  cross  profiles  at  the  high-bank  positions.  Photo 
interpreters  interpolated  elevations  along  the  high  banks  between  those 
data  points  with  the  aerial  photos  as  an  aid.  In  attempt  to  use  avail- 
able contour  maps  was  unsuccessful,  since  their  data,  when  available, 
were  much  coarser  than  those  which  the  photo  interpreters  could 
reasonably  interpret.  Most  side  channels  exhibited  a smooth  change  in 


elevation  along  the  high  banks , typically  changing  a few  feet  per  mile, 
thus  making  interpolation  simple.  Since  the  parameter  values  were 
calculated  as  a function  of  water  elevation  up  to  a maximum  elevation 
equal  to  that  at  which  50  percent  of  the  high  banks  were  overtopped 
(see  Part  III),  the  results  of  this  study  are  independent  of  elevations 
along  a high-bank  bluff. 

28.  The  water  elevations  of  standing  pools  were  interpolated  by 
studying  the  shapes  of  the  pools  and  attempting  to  locate  the  inter- 
section of  the  water  surface  with  cross  profiles.  The  distances  of  the 
land-water  interfaces  from  the  high  banks,  as  well  as  the  pool  width, 
were  measured  from  the  aerial  photography  and  scaled  onto  the  cross 
profiles  to  provide  the  approximate  intersection  points.  Locating  a 
pool's  elevation  was  not  difficult  in  most  cases,  since  the  side-channel 
shapes  did  not  change  radically  in  the  time  between  aerial  photo 
coverage  and  measurement  of  profiles. 

Topography  Simulation 

29.  The  production  of  topographic  data  that  describe  the 
geometry  of  each  side  channel  required  five  operations  and  resulted  in 
three  products.  The  five  operations  were  the  manual  digitizing  of  the 
data  on  the  base  maps  and  profile  plots  in  a format  suitable  for  computer 
operations,  production  of  a set  of  (XYZ)  data  from  the  base  map  and 
profile  data,  computer  calculation  of  an  elevation  grid  array  over  the 
side  channel,  computer  calculation  of  a contour  map  of  the  side  channel, 
and  the  manual  digitizing  of  that  contour  map  to  produce  a file  contain- 
ing the  (XYZ)  coordinates  of  the  contour  lines. 

30.  The  three  products  resulting  from  this  sequence  of  operations 
were  the  elevation  grid  array,  the  contour  map,  and  the  computer  file  of 
(XYZ)  coordinates  of  the  contours.  The  elevation  grid  array  and  the 
computer  file  of  (XYZ)  coordinates  were  produced  to  provide  topographic 
data  in  a form  necessary  for  calculation  of  the  parameter  values.  The 
contour  map  was  an  intermediate  product  in  producing  the  contour  (XYZ) 
coordinates . 


17 


r 


Data  retrieval 

31.  The  data  on  the  base  map  and  fathometer-survey  plots  were 
digitized  using  the  same  equipment  described  in  paragraph  15.  The  infor- 
mation retrieved  from  the  base  map  included  the  (XY)  coordinates  for  the 
boundary  of  the  side  channel,  the  base  map  scale,  the  identification 
numbers  and  (XY)  coordinates  for  the  end  points  of  each  cross  profile, 
the  identification  numbers  and  (XY)  coordinates  for  each  section  of  the 
thalweg  profile  where  it  intersected  the  cross  profiles,  and  the  (XYZ) 
coordinates  of  all  elevation  data,  such  as  waterlines,  high  banks,  and 
dikes.  The  information  retrieved  from  each  of  the  cross-profile  plots 
included  the  elevation  scale,  the  mean  sea  level  elevation  of  the  water 
surface,  the  identification  number,  the  high-bank  positions , the  orienta- 
tion of  the  profile  relative  to  direction  of  water  flow,  and  the  (XY) 
coordinates  of  the  profile  in  sufficient  detail  to  describe  the  profile 
with  straight  lines  between  data  points.  The  information  retrieved 

from  each  thalweg  profile  included  the  elevation  scale , the  mean  sea 
level  elevation  of  the  water  surface,  the  (XY)  coordinates  of  the  profile 
in  sufficient  detail  to  describe  the  profile  with  straight  lines  between 
data  points  , and  the  location  of  each  position  on  the  profile  where  it 
intersected  a cross  profile. 

32.  Because  of  the  differences  in  side-channel  sizes  and  com- 
plexity of  data,  digitizing  time  varied  from  2 to  6 hr/per  side  channel. 
Typically,  a relatively  unskilled  equipment  operator  could  digitize  all 
the  necessary  data  for  one  side  channel  in  an  average  of  about  4 hr. 

At  the  conclusion  of  the  digitizing  process,  the  digitized  data  were 
on  magnetic  tape  and  ready  for  use  in  the  next  operation. 

(XYZ)  data  preparation 

33.  The  data  for  a side  channel  contained  on  a magnetic  tape  from 
the  digitizer  were  used  as  input  to  a computer  program  that  produced  a 
second  data  file  containing  (XY)  coordinates  of  the  boundary  (high  banks 
and  ends)  of  the  side  channel  and  (XYZ)  coordinates  of  the  elevations 
within  the  side  channel. 

34.  The  coordinates  of  the  boundary  and  the  (XYZ)  coordinates  of 
the  elevation  data  derived  from  the  base  map  were  simply  repeated  in  the 
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output  data  file  by  that  program.  The  program  then  "recognized"  each 
profile  by  its  identification  number  and  superimposed  it  on  the  other 
base  map  data  in  the  form  of  (XYZ)  coordinates  between  the  positions  of 
the  high-bank  coordinates  for  that  profile.  Each  profile  was  positioned 
properly  across  the  side  channel  relative  to  the  direction  of  water 
flow  and  scaled  linearly  so  that  the  high-bank  positions  on  the  profile 
were  made  to  coincide  with  the  high-bank  end  positions  for  that  profile 
on  the  base  map.  Finally,  each  digitized  point  along  the  profile  was 
transformed  into  an  (XYZ)  coordinate  on  the  base  map. 

35-  The  same  operation  was  performed  with  the  sections  of  the 
thalweg  profile,  except  the  sections  were  scaled  linearly  so  that  the 
end  points  of  each  section  were  made  to  coincide  with  the  points  of 
intersection  between  the  thalweg  and  the  cross  profiles. 

Grid  process 

3 6.  The  set  of  (XY)  points  defining  a side-channel  boundary  and 
the  set  of  (XYZ)  coordinates  within  that  boundary  were  input  to  a com- 
puter program  that  calculated  elevations  on  a grid  across  that  side 
channel.  "Calculated  elevations  on  a grid"  means  that  an  elevation  was 
calculated  at  each  intersection  of  equally  spaced  horizontal  and  vertical 
lines  (20-m  spacing)  superimposed  over  the  side  channel.  The  (XY) 
coordinates  defining  the  side-channel  boundary  were  used  by  the  program 
to  restrict  its  calculations  to  the  region  within  the  side-channel  high 
banks  and  ends . 

37-  The  elevation  grid  was  calculated  because  data  in  that  format 
make  many  analysis  procedures  straightforward,  including  calculations  of 
values  of  some  of  the  parameters  of  this  study  and  the  calculation  of  the 
contour  maps.  (The  calculations  of  parameter  values  and  the  contour  maps 
are  discussed  in  Part  III.)  Also,  the  maximum  amount  of  informational 
content  can  be  compressed  into  the  computer  data  file  compared  with  the 
space  required  by  other  (than  grid)  formats,  since  only  the  uniform 
grid  spacing  and  the  elevations  of  the  grid  points  are  required  to 
describe  the  geometry  rather  than  a set  of  (XYZ)  coordinates.  The  eleva- 
tion at  any  (XY)  position  can  be  found  by  retrieving  the  elevation 
from  the  proper  row  and  column  position. 


38.  The  procedure  used  for  calculating  the  elevation  at  each  grid 
position  is  one  that  has  been  sucessfully  applied  at  the  WES  in  many 
prior  studies  involving  reconstruction  of  topographic  surfaces  from 
sparse  data.  The  calculational  procedure  is  as  follows:  A local  coordi- 

nate system  is  first  centered  at  the  grid  position  for  which  an  eleva- 
tion is  to  be  calculated,  and  the  space  about  the  grid  position  is 
divided  into  quadrants.  The  data  point  in  each  quadrant  closest  to  the 
grid  position  (nearest  neighbor)  is  selected  out  of  the  total  available 
data  set.  The  elevations  and  positions  of  those  four  data  points  are 
used  to  calculate  the  elevation  at  the  grid  position,  using  the  inverse- 
distance-  square  weighted  elevations  of  the  four  data  points.  The 
algortihm  used  is 


where 

Z = the  elevation  at  the  grid  position 

i = the  index  over  the  four  nearest  neighbors 
= the  elevation  of  the  ith  nearest  neighbor 

R.  = the  distance  of  the  i^*1  nearest  neighbor  from  the  grid 
position 

If  a nearest  neighbor  cannot  be  found  in  all  quadrants,  as  frequently 
occurs  when  a grid  position  is  close  to  the  boundary  of  the  site,  as 
many  nearest  neighbors  as  are  found  are  used.  If  a data  point  is 
located  within  the  immediate  vicinity  of  the  grid  position  (within  a 
radius  less  then  one-tenth  the  space  between  grid  positions),  the  eleva- 
tion of  that  data  point  is  used  for  the  grid  position  elevation. 

39-  A topographic  surface  produced  with  the  algorithm  above  has 
the  qualities  of  being  smoothly  varying  with  no  discontinuties  and  of 
providing  an  exact  fit  in  the  locale  of  all  data  points . Other 
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advantages  of  the  grid  representation  are  that  the  interpolative 
procedures  bring  out  surface  details  not  immediately  apparent  and 
frequently  overlooked  when  the  data  are  handled  by  other  procedures . 

40.  The  informational  content  of  a grid  array  representation  of  a 
topographic  surface  is  a function  of  the  grid  spacing.  Since  the  grid 
array  representation  is  achieved  by  an  interpolative  procedure,  it  can 
contain  no  more  information  about  the  surface  than  the  set  of  (XYZ)  coor- 
dinates used  as  input.  Based  on  past  experience,  a grid  spacing  of  20  m 
was  chosen  for  the  grid  arrays  of  all  side  channels  based  on  the  esti- 
mated accuracy  and  distribution  of  the  data,  the  resolution  required  for 
the  calculation  of  parameter  values,  the  desire  to  make  any  errors  in  the 
final  results  due  to  the  grid  spacing  both  systematic  and  in  the  same 
direction  for  results  for  all  side  channels , and  computer  time  and 
memory  space  limitations. 

Contour  maps 

1*1.  Computer-calculated  and  -plotted  contour  maps  were  produced 
for  each  of  the  18  side  channels.  Reproductions  of  these  maps  are  con- 
tained in  Volume  II.  The  contour  maps  were  produced  to  yield  a simple 
and  rapid  means  of  checking  the  contents  of  the  grid  array  file  for 
errors  in  the  data  input  to  the  grid  array  calculation,  and  to  provide 
the  topographic  data  in  a form  needed  for  a subsequent  operation.  Some 
input  data  errors  were  discovered  using  this  process  and  subsequently 
corrected. 

42.  The  computer  procedure  used  for  calculating  the  profiles  was 
as  follows:  The  topographic  surface  of  a side  channel  was  broken  into 

a series  of  grid  squares,  where  each  grid  square  was  defined  by  the 
elevations  of  the  grid  positions  at  its  four  corners.  Each  grid  square 
was  further  subdivided  by  cutting  it  diagonally  into  two  triangles,  and 
then  planes  were  uniquely  fit  to  the  two  sets  of  triangularly  arranged 
points  in  each  grid  square.  A series  of  equally  separated  horizontal 
planes  was  constructed  starting  at  mean  sea  level  with  a separation 
equal  to  the  contour  interval  (5  ft).  Whenever  a "triangular" 
topographic  surface  plane  intersected  a horizontal  plane,  the  coordinates 
of  the  line  of  intersection  were  calculated  and  computer  plotted.  Since 


21 


the  surface  representation  as  a grid  array  was  smoothly  varying  with  no 
discontinuities,  the  total  series  of  lines  of  intersection  calculated 
and  computer  plotted  formed  closed  contour  lines  on  completion  of  the 
process  over  the  entire  grid  array. 

U3.  Since  the  contour  lines  were  calculated  as  intersections  of 
planes  with  planes,  each  contour  line  was  a series  of  short,  straight- 
line  segments  as  seen  in  the  maps  in  Volume  II.  The  contour  maps  were 
plotted  with  straight-line  segments  between  calculated  positions  on 
the  contour  lines  rather  than  with  smoothed  lines,  because  smoothing 
techniques  ordinarily  introduce  errors  into  the  contour  map  in  the 
form  of  contour  line  dislocation. 

kk.  A series  of  contour  maps  was  plotted  for  selected  side 
channels  at  several  contour  intervals  to  determine  the  smallest  interval 
that  was  consistent  with  the  data  content  of  the  grid  array.  It  was 
originally  intended  to  produce  maps  with  2-ft  contour  intervals;  however, 
such  maps  showed  distinctive  effects  of  interpolating  between  larger 
intervals,  i.e.  there  was  a high  occurrence  of  parallel  contour  lines 
over  most  substructures  of  a bottom  surface.  Therefore,  5-ft  intervals 
were  chosen,  since  maps  plotted  at  that  and  larger  intervals  did  not 
show  the  above-mentioned  effect  to  any  appreciable  degree. 

Contour  (XYZ)  coordinates 

1+5.  Contour  data  calculated  from  the  grid  array  were  used  to  plot 
contour  maps  automatically.  The  maps  (shown  in  Volume  II)  were  then 
digitized  manually  with  the  same  equipment  previously  described  to  pro- 
vide the  topographic  data  in  a form  suitable  for  calculating  parameters 
d_  and  e_,  paragraph  5*  While  manual  digitizing  was  laborious  and  time- 
consuming,  the  total  project  time  and  funds  did  not  permit  the  develop- 
ment of  a more  versatile  automated  procedure.*  The  manual  procedure  for 
all  maps  was  the  same  and  consisted  of  digitizing  each  contour  line  as 
a closed  loop,  thereby  providing  the  elevation  and  coordinates  of  each 
line  in  a computer  file  for  use  in  calculating  some  of  the  parameters. 

* A new  procedure  for  performing  this  step  with  increased  accuracy  and 
a drastic  reduction  in  time  and  cost  per  side  channel  was  developed 
shortly  after  project  completion. 
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PART  III:  CALCULATIONS  AND  RESULTS 
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1*6.  The  parameters  for  which  values  were  calculated  are  listed 

in  paragraph  5,  hut  are  repeated  below  for  convenience. 

a_.  Center-line  length. 

b^  Average  width  between  high  banks. 

c_.  Water  volume  as  a function  of  water  elevation. 

d_.  Shoreline  length  as  a function  of  water  elevation. 

e_.  Water  surface  area  as  a function  of  water  elevation. 

_f.  Shoreline  development  as  a function  of  water  elevation. 

Rate  of  change  of  water  surface  area  with  respect  to  water 
elevation  (derivative  of  water  surface  area  with  respect 
to  water  elevation)  as  a function  of  water  elevation. 

h.  Ratio  of  water  surface  area  to  volume  as  a function  of 
water  elevation. 

i_.  Ratio  of  shoreline  length  to  water  surface  area  as  a 
function  of  water  elevation. 

j_.  Bottom  surface  area  underwater  as  a function  of  water 
elevation  and  water  depth. 

k.  Water  cross-sectional  area  as  a function  of  water  eleva- 
tion at  selected  sampling  locations  (stations). 


Water  Elevation 


1*7.  "Water  elevation"  is  defined  as  the  elevation  of  the  water 
within  the  side  channel.  It  is  important  to  recognize  certain  problems 
associated  with  specifying  water  elevation,  the  errors  that  can  arise  in 
the  calculational  results  because  of  those  problems,  and  the  error 
correction  procedures  built  into  the  procedures  for  evaluating  the 
parameters . 

1*8.  Three  conditions  associated  with  water  surface  slope  and 
impoundment  affect  the  meaning  and  assignment  of  water  elevation: 

a.  There  is  a difference  between  the  water  elevation  in  the 
side  channel  and  that  in  the  main  channel. 

b.  There  is  a change  in  water  elevation  in  the  side  channel 
relative  to  mean  sea  level  due  to  the  surface  slope  of 
moving  water  in  the  side  channel. 


c_.  Any  impounded  water  within  a side  channel  is  typically  not 
at  the  same  elevation  as  that  in  the  main  channel  or  other 
parts  of  the  side  channel. 

Each  of  these  conditions  is  discussed  below. 

U9 - Since  only  the  main  channel  water  elevation  was  available  at 
the  time  the  fathometer  profiles  were  measured  and  at  the  time  of  aerial 
photo  coverage,  and  no  data  were  available  to  specify  the  difference 
between  main-  and  side-channel  water  elevations,  the  elevation  of  the 
latter  was  accepted  as  equal  to  the  former  when  they  were  connected. 

In  addition,  the  difference  in  water  levels  is  not  constant,  but  varies 
uniquely  in  each  side  channel  with  river  stage,  so  that  the  difference 
at  the  time  of  photo  coverage  was  not  the  same  as  that  at  the  time  of 
fathometer  measurements,  nor  the  same  from  one  side  channel  to  the  next.- 
This  problem  does  not  affect  the  shape  of  the  parameters  (how  the  values 
of  the  parameters,  which  are  a function  of  water  elevation,  change  with 
changes  in  water  elevation),  but  rather  introduces  errors  into  the 
absolute  values  of  those  parameters.  Since  the  difference  in  water 
elevation  referred  to  is  considered  small  (typically  less  than  1 ft  for 
any  side  channel  at  any  time),  the  introduced  error  is  considered  negli- 
gible compared  with  errors  due  to  sparseness  of  data. 

50.  Because  of  the  slope  of  the  water  surface  in  a side  channel, 
the  water  elevation  relative  to  a fixed  datum  varies  along  the  side 
channel  when  water  is  flowing  through  it,  but  is  almost  constant  when 
flow  is  stopped.  The  fathometer  profiles  were  measured  at  a time  of 
maximum  flow — when  the  water  surface  slope  within  the  side  channel  most 
closely  approximated  that  within  the  main  channel.  The  fathometer 
profiles  contain  the  elevations  affected  by  water  surface  slope  at  that 
time,  since  all  elevations  in  those  data  are  relative  to  the  water 
surface.  The  values  of  the  parameters  are,  therefore,  already  partially 
adjusted  for  that  effect  for  the  water  elevations  when  water  is  flowing 
through  the  side  channels.  The  end-to-end  difference  in  water  elevation 
is  typically  no  more  than  1.5  ft  for  any  side  channel,  which  is  small 
compared  with  errors  introduced  in  interpretation  of  the  fathometer  and 
aerial  photo  data. 
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51.  All  parameters  in  paragraph  46,  except  ji  and  b,  were  calcu- 
lated as  a function  of  water  elevation  at  a single  location,  recognizing 
the  effects  of  water  surface  slope  as  noted  above.  That  location  was 
the  approximate  lengthwise  center  of  the  side  channel.  The  water  eleva- 
tion in  the  main  channel  opposite  that  location  was  found  (rounded  off  to 
the  nearest  foot)  from  interpolating  between  the  mean  sea  level  eleva- 
tions available  on  the  aerial  photography  and  the  elevation  assigned  to 
water  in  the  side  channel  connected  with  the  main  channel . The  results 
presented  in  this  report  as  a function  of  water  elevation  are  relative 

to  that  center  point.  The  single  exception  is  parameter  k,  for  which 
values  were  referenced  to  the  water  elevation  at  the  location  of  the 
cross-sectional  calculation  (see  paragraphs  76  and  77). 

52.  As  the  water  level  drops  in  a side  channel,  various  struc- 
tures, particularly  sandbars  at  the  ends  of  the  side  channels  and  dikes 
within  or  at  the  side-channel  ends,  restrict  water  flow.  As  the  water 
level  continues  to  drop,  pools  are  formed  with  different  water  eleva- 
tions. The  water  level  within  a pool  does  not  normally  remain  constant, 
however.  The  local  climatic  conditions,  the  hydraulic  head,  the 
composition  of  the  dike  or  other  control  structures,  and  the  composition 
of  the  base  material  affect  the  percolation,  evaporation,  and  seepage 
rates  in  lowering  the  pool  level.  Local  rainfall,  runoff,  and  fluctu- 
ations of  the  main  channel  above  the  elevations  of  control  structures 
work  to  raise  the  pool  level.  Finally,  a time  comes  when  the  eleva- 
tion of  the  main  river  rises  above  the  elevation  of  one  of  the  control 
structures,  and  the  pool  is  again  connected  to  the  main  channel.  The 
locations  and  elevations  of  the  major  control  structures  are  known,  but 
there  is  no  general  theoretical  relation  developed  that  takes  all 
environmental  controlling  conditions  into  account,  nor  are  empirical  data 
generally  available  to  describe  the  elevations  of  the  impounded  water 
relative  to  the  elevation  of  the  main  river  channel. 

53-  If  that  information  were  available,  each  of  the  parameters 
evaluated  in  this  study  would  be  double-valued  over  the  water  elevation 
range,  from  that  elevation  at  which  water  was  impounded  anywhere  in  the 
side  channel  down  to  the  minimum  water  elevation  in  that  side  channel. 
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Lacking  those  data,  the  parameters  were  calculated  for  the  entire  side 
channel,  ignoring  the  impoundment  effects.  The  values  of  the  parameters 
calculated  in  this  study  for  water  elevations  lower  than  an  " impoundment 
elevation"  for  a side  channel  are  between  the  double  values.  There  is 
no  assurance  that  a single-valued  function  is  at  the  mean  of  the  double- 
valued function  for  any  of  the  parameters. 

54 . The  state-of-the-art  incapability  to  provide  a means  of  cal- 
culating the  double-  rather  than  single-valued  functions  over  the  range 
of  impoundment  is  not  particularly  distressing,  considering  the  purpose 
of  this  study.  The  purpose  is  to  provide  geometric  characteristics  (as 
one  part  of  the  physical  characteristics)  for  correlation  with  wildlife 
population  structures.  The  accuracy  of  such  a correlation  depends  on  the 
appropriateness  of  the  parameters  chosen  to  quantitatively  describe  the 
characteristics  and  the  accuracy  with  which  the  values  of  those  param- 
eters are  calculated.  It  is  well  known  mathematically  that  the  accuracy 
of  the  results  of  a correlation  is  a function  of  the  accuracy  of  the 
parameters  used  in  the  correlation,  and  that  the  correlation  accuracy 
cannot  be  greater  than  the  worst  accuracy  of  any  of  the  set  of  parameters 
used  in  the  correlation.  The  accuracy  of  the  parameters  calculated  in 

this  study  is  much  greater  than  that  of  other  than  geometric  parameters  i 

or  of  the  biological  data  and  is,  therefore,  more  than  sufficient  for 

the  intended  purpose.  The  reason  for  the  difference  in  accuracy  is  5 

primarily  one  of  sample  size.  Expressed  simply,  the  geometric  parameter 
values  are  typically  based  on  several  thousand  "samples"  per  side 

' 

channel,  while  such  sample  sizes  are  not  normally  obtainable,  for 
example,  for  animal  populations. 

Pool  Parameters 

55.  In  addition  to  calculation  of  the  parameters  over  the  entire 
range  of  water  elevations  of  all  side  channels,  side  channels  Osborne, 

Fort  Chartres,  and  Liberty  were  selected  for  more  detailed  calculations.  | 

The  levels  and  locations  of  the  water-control  structures  were  identified, 
and  all  parameters  were  evaluated  for  each  separate  "pool"  formed  by  the 


26 


control  structures  as  the  water  elevation  decreased.  Pool  1 of  all  three 


side  channels  consisted  of  the  entire  side  channel  for  the  range  of  water 
elevations  between  the  maximum  (see  paragraph  5 6)  and  the  elevation  at 
which  any  impoundment  caused  by  a control  structure  first  took  place, 
forming  pools  2 and  3.  When  a second  control  structure  was  located  in 
pools  1 or  2,  the  parameters  were  calculated  between  the  upper  water 
elevation  level  at  which  the  pool  was  formed  and  the  lower  elevation  at 
which  it  separated  into  pools  4 and  5.  When  a pool  did  not  contain  any 
additional  control  structures  as  the  water  elevation  decreased,  the 
parameters  were  calculated  between  the  upper  water  elevation  level  at 
which  the  pool  was  formed  and  the  water  elevation  at  which  negligible 
water  remained  in  the  pool.  The  three  selected  side  channels  con- 
tained 5,  3,  and  5 pools,  respectively.  The  locations  of  the  pools 
within  the  side  channels  are  shown  in  fig.  6.  Fig.  6 can  be  understood 
by  viewing  it  in  conjunction  with  the  contour  maps  for  the  above- 
mentioned  side  channels  in  Volume  II.  The  dikes,  which  were  the  con- 
trol structures  for  those  side  channels,  are  also  shown  on  the  contour 
maps  in  Volume  II.  A review  of  fig.  2 would  also  be  helpful,  since  it 
shows  Osborne  at  low  water,  revealing  the  dikes  that  formed  the  pools. 
Range  of  calculations 

5 6.  Calculations  were  performed  at  2-ft  water  elevation  intervals 
for  those  parameters  that  were  functions  of  water  elevation,  starting 

at  the  water  elevation  at  which  the  side  channel  was  dry  or  nearly  dry 
and  continuing  up  to  the  elevation  at  which  50  percent  of  the  high  banks 
was  overtopped. 

57.  Parameters  c_,  and  k_  (paragraph  46)  were  calculated  directly 
from  the  grid  array  data  at  each  2-ft  interval . Parameters  d_  and  e.  were 
calculated  using  the  digitized  contour  data,  which  were  at  5-ft 
intervals.  The  2-ft-interval  data  for  the  parameters  were  interpolated 
from  the  5-ft  data.  The  total  calculational  procedure  for  all  parameters 
is  described  in  paragraphs  61-77. 

58.  The  upper  and  lower  water-level  extremes  used  in  the  calcula- 
tions for  a specific  side  channel  were  chosen  by  inspection  of  the  base 
map  and  the  contour  map,  respectively.  Since  the  side-channel  high  banks 


b.  Three  pools  of  side  channel  Fort  Chartres 

Fig.  6.  Locations  of  pools  in  selected  side 
channels  (sheet  1 of  2) 


Five  pools  of  side  channel  Liberty 
Fig.  6.  (sheet  2 of  2) 


were  typically  smoothly  varying,  it  was  possible  to  choose  the  upper 
water-level  extreme  easily  and  unambiguously.  Selection  of  the  lower 
extreme  was  more  subjective.  The  general  rule  was  that  the  extreme  be 
picked  at  that  elevation  at  which  separated  standing  pools  of  no  signif- 
icant size  (compared  with  the  to+al  side-channel  extent)  remained. 

59-  When  calculations  were  performed  for  those  levels  at  which  a 
portion  of  the  high  bank  was  overtopped,  a vertical  wall  was  assumed  at 
the  locations  of  the  overflow,  so  that  the  parameters  could  be  evaluated 
up  to  the  elevation  of  50  percent  overtop. 

Water-elevation  reference 

60.  The  water  levels  at  which  the  parameter  values  were  calculated 
were  the  mean  sea  levels  of  the  water  at  the  center  of  each  side  channel 
(see  paragraph  51).  For  ease  in  interpretation,  all  water  elevations 
reported  in  the  results  were  also  converted  to  the  St.  Louis  (Market 
Street)  equivalent  gage.* 


Calculational  Procedures 


61 . The  procedures  for  computer  calculation  of  each  parameter 
listed  in  paragraph  U6  were  identical  for  all  side  channels  and  are 
described  below. 

Parameter  a:  Center-line  length 

62.  The  center-line  length  was  calculated  from  the  base  map.  A 
center  line  midway  between  high  banks  of  the  side  channel  was  placed  on 
the  photo-interpreted  representation  of  the  side  channel  on  the  base  map. 


* The  water  elevation  at  any  position  along  the  Mississippi  River  is 
normally  expressed  relative  to  the  scale  of  the  gage  closest  to  that 
position,  or  relative  to  mean  sea  level.  It  is  difficult  for  a person 
to  interpret  the  significance  of  a gage  reading  unless  he  has  experi- 
ence with  that  gage.  A somewhat  artificial  gage  system,  the  St.  Louis 
(Market  Street)  equivalent  gage,  was  devised  and  is  used  by  the 
St.  Louis  District  to  help  District  personnel  better  relate  to  river 
stage  readings  at  locations  where  they  do  not  have  experience.  Using 
that  system,  the  water  elevation  at  any  location  on  the  river  between 
St.  Louis,  Missouri,  and  Cairo,  Illinois,  is  converted  to  the  equiv- 
alent elevation  that  would  be  read  on  the  St.  Louis  Market  Street 
gage  if  the  river  was  in  a steady-state  condition. 


Length  was  measured  by  digitizing  the  center  line  as  a series  of  very 
short,  straight  segments,  and  calculating  the  total  center-line  length 
as  the  sum  of  the  segment  lengths  with  the  following  equation: 


X.)2  + (Y  - Y.)2 
x 1+1  1 


L = the  center-line  length 
N = the  total  number  of  segments 

X.,Y.  = the  Cartesian  coordinates  of  the  i^  digitizer  point  on  the 
center  line 

(N+l)  = the  number  of  digitized  points  defining  the  center  line 


The  straight-line  segments  were  made  small  enough  so  that  the  total 
digitized  center  line  was  essentially  continuous. 

63.  When  a side  channel  had  branches  (see  contour  map  II-5  in 
Volume  II)  so  that  more  than  one  center  line  coulu  be  defined,  the 
length  of  the  longest  center  line  was  calculated. 

oU . The  "ends"  of  a side  channel  on  the  base  map  were  not  neces- 
sarily at  the  juncture  of  that  side  channel  with  the  main  channel,  but 
were  generally  fixed  at  the  first  structure  that  controlled  water  flow 
within  the  side  channel  from  that  juncture. 

Parameter  b : Average 

width  between  high  banks 

65.  The  average  width  between  high  banks  was  calculated  from  the 
base  map.  On  the  average,  10  lines  were  constructed  perpendicular  to 
the  center  line,  and  the  distances  between  high  banks  on  those  lines 
were  calculated  and  averaged  with  the  following  equation: 


it 

I V(x2i  - xu>2  + (Ia  - V1 
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where 


W = the  average  side-channel  width 
N = the  total  number  of  lines 


the  coordinates  of  one  high-bank  position  on  the  i line 

the  coordinates  of  the  other  high-bank  position  on  the  i '' 
line 


Parameter  c : Water  volume 

as  a function  of  water  elevation 


66.  The  water  volume  as  a function  of  water  elevation  was  calcu- 
lated using  the  elevation  grid  array  as  a data  source.  The  volume  for  a 
specific  water  elevation  was  calculated  by  summing  the  volumes  of  water 
over  all  grid  squares,  where  a grid  square  was  that  section  of  the  side- 
channel  bottom  surface  determined  by  four  adjacent  grid  elevations.  The 
volume  of  water  over  any  grid  square  was  calculated  in  the  following 
manner:  The  grid  square  was  first  cut  into  two  triangles  with  a vertical 

plane  passing  through  a grid  square  diagonally.  The  vertical  plane  also 
cut  the  rectangular  column  of  water  above  that  grid  square  into  two 
triangular  columns.  The  volumes  of  both  triangular  columns  bounded  by 
the  plane  of  the  water  surface,  the  plane  formed  by  the  triangle  on  the 
bottom,  and  vertical  planes  passing  through  the  three  sets  of  vertices 
of  the  triangle  were  calculated  and  summed  to  yield  the  volume  of  water 
over  that  grid  square.  A mathematical  expression  for  the  entire  opera- 
tion is  as  follows: 


V = the  volume  of  water  in  the  side  channel  when  water  elevation 
is  at  the  k*'*1  level 

M = the  number  of  grid  squares  in  the  Y direction 
N = the  number  of  grid  squares  in  the  X direction 
D = the  horizontal  distance  between  elevation  grid  points 


= the  elevation  of  water  at  the  k 


level 


— H ti»h 

ZV . = the  average  of  the  three  grid  elevations  in  the  R triangle 

* J 

3 

of  the  (i,j)  grid  square  = V Z..  /3 

p=L  1JP 

Z1? . = elevation  of  the  grid  point  in  the  R*-'1  triangle  of 

1J p the  (i,j)th  grid  square 

and  where  the  summation  is  completed  only  for  those  triangles  of  bottom 

surface  within  the  grid  squares  where  the  water  elevation  is  above  the 

— R 

bottom  surface  (E  > Z ). 

K 1 J 

Parameter  d:  Shoreline  length 

as  a function  of  water  elevation 


67.  The  shoreline  length  as  a function  of  water  elevation  was 
calculated  from  the  digitized  contour-line  data.  The  shoreline  length 
at  a specific  water  elevation  was  calculated  by  summing  the  lengths  of 
all  contour  lines  (as  closed  loops)  at  that  elevation.  Since  the 
contour  lines  were  at  5-ft  intervals,  a table  was  formed  of  shoreline 
length  at  5-ft  water-elevation  intervals.  The  mathematical  procedure 
used  was  as  follows : 


V'L  Vi 

-11 

j=l  i=l 


<ck 

i+l.J 


i+l.J 


.k  \ 

i j / ’ 


= xk 

N+l,j  lj 


Yk  = Xk 
N+l.J  1J 


where 


L = the  shoreline  length  when  water  is  at  the  k level 

k th 

M = the  number  of  contour  lines  with  water  at  the  k“  level 

N = the  number  of  digitized  data  points  on  the  j11*  contour 
line  with  an  elevation  at  the  k^h  level 

(xk.,Yk.j  = the  coordinates  of  the  i*'11  point  on  the  eontov 

1J  that  has  an  elevation  at  the  k^h  level 


contour  line 


and  where  the  contour  lines  were  "closed"  in  the  calculation  by  connect- 


ing the  coordinates  of  the  first  point  in  the  string  of  numbers  with  the 
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coordinates  of  the  last  point  by  specifying  the  condition 


(4*1, ij-  (4r4j) 


Parameter  e:  Water  surface  area 

as  a function  of  water  elevation 

68.  The  water  surface  area  as  a function  of  water  elevation  was 
calculated  from  the  digitized  contour-line  data.  The  surface  area  at 
any  water  elevation  was  calculated  by  summing  the  areas  inclosed  by  all 
contour  lines  with  that  elevation  which  were  covered  with  water.  The 
area  inclosed  by  a contour  line  was  calculated  using  the  trapezoidal 
rule  on  the  set  of  coordinates  defining  that  contour  line  as  a closed 
loop.  The  mathematical  procedure  was  as  follows: 


\ - I 1 2 (4.i,j  - 4j)(4*i,j  * 4j) 


J=1  i=l 


X* 

N+1,J 


N+l.J 


where 


A,  = the  water  surface  area  when  water  is  at  the  k 
k . . 

th 

M = the  number  of  contour  lines  at  the  k level 


(4r4j) 


the  number  of  digitized  data  points  on  the 
line  with  an  elevation  at  the  k'*'*1  level 

the  coordinates  of  the  i^  point  on  the  j1 
that  has  an  elevation  at  the  k^h  level 


1 evel 


contour 


contour  line 


and  where  the  condition  , . ,Y^  , .)  = [X.  satisfies  the  condi- 

V N+l , j N+l,j/  V lj’  lj/ 


tion  of  contour-line  closure. 

Parameter  f:  Shoreline  develop- 

ment as  a function  of  water  elevation 


69.  Shoreline  development  is  a parameter  commonly  used  as  a 
measure  of  a water  body's  circularity.  It  can  be  easily  shown  that  the 
minimum  extremum  circumference  of  a two-dimensional  geometric  figure 
inclosing  a region  must  have  a circular  shape.  The  value  of  shoreline 
development  of  a circular  water  body  is  unity,  whereas  that  of  all  water 


bodies  with  other  geometric  configurations  is  greater  than  unity.  Thus, 
the  departure  of  a water  body's  shoreline  development  value  from  unity 
is  a measure  of  its  "differentness"  from  a circle.  The  parameter  values 
were  calculated  using  the  following  equation: 


D 


k 


where 

th 

D = the  shoreline  development  when  water  is  at  the  k level 
K 

L = the  shoreline  length  when  water  is  at  the  k level 

^ th 

A^  = the  water  surface  area  when  water  is  at  the  k level 


Parameter  g:  Rate  of 

change  of  water  surface  area 
with  respect  to  water  elevation 

TO.  The  rate  of  change  (derivative)  of  water  surface  area  with 
respect  to  water  elevation  was  calculated  using  the  previously  calculated 
surface  area  data.  A smooth  curve,  a third-order  polynomial,*  was 
analytically  fit  to  the  surface  area  as  a function  of  water  elevation 
data  in  the  immediate  locale  of  each  water  elevation  for  which  a deriva- 
tive would  be  calculated.  First  derivatives  were  then  calculated  at 
the  2-ft-interval  water  elevations  from  the  analytic  expressions  for  the 
data  (one  analytic  expression  per  data  point). 

Parameter  h:  Ratio  of  water 

surface  area  to  volume  as 
a function  of  water  elevation 

71.  The  ratio  of  water  surface  area  to  volume  as  a function  of 
water  elevation  was  calculated  using  the  previously  calculated  shoreline- 
length  and  water-surface-area  data.  Ratios  were  taken  at  each  water 
elevation  as  expressed  by  the  following  relation: 


* A second-order  polynomial  was  used  for  data  points  close  to  the  low- 
water  and  high-water  ends  of  the  data. 
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where 


the  ratio  of  water  surface  area  to  water  volume  when  water 
is  at  the  kth  level 

Vi 

the  water  surface  area  when  the  water  is  at  the  k n level 
the  water  volume  when  water  is  at  the  k^  level 

Ratio  of  shore- 
line length  to  water  surface 
area  as  a function  of  water  elevation 

72.  The  ratio  of  shoreline  length  to  water  surface  area  as  a 
function  of  water  elevation  was  calculated  using  the  previously  calcu- 
lated shoreline-length  and  water -surf ace -are a data.  Ratios  were  taken 
at  each  water  elevation  as  expressed  by  the  following  equation: 


where 

R^  = the  ratio  of  water  surface  area  to  water  volume  when  the 
water  is  at  the  kth  level 

k "V  Vi 

L = the  shoreline  length  when  the  water  is  at  the  k level 
k -v  Vi 

A = the  water  surface  area  when  the  water  is  at  the  k level 

Parameter  j : Bottom  surface 

area  underwater  as  a function 
of  water  elevation  and  water  depth 

73.  The  relation  between  water  depth  and  wildlife  populations  is 
not  completely  understood.  There  is  information  (mostly  qualitative), 
however,  regarding  the  "typical"  water  depths  acceptable  to  different 
waterfowl  and  fish  species,  and  acceptable  to  fish  species  during  the 
major  phases  of  their  life  cycles.  Using  that  information,  seven  depth 
class  ranges  were  defined  as  follows, 
si.  Surface  to  2 ft. 
b . From  2 to  U f t . 


Parameter  i : 
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I 


J 


1 


c_.  From  h to  6 ft. 
d . From  6 to  10  ft . 
e_.  From  10  to  15  ft. 
f_.  From  15  to  20  ft . 

£.  In  excess  of  20  ft. 

An  attempt  was  made  to  class  the  depths  such  that  conditions  that  influ- 
ence wildlife  population  densities  (for  this  parameter)  are  uniform  in 
each  depth  class  range,  hut  differ  from  one  class  to  another. 

74.  The  data  for  the  bottom  surface  area  under  water  as  a function 
of  water  elevation  and  water  depth  were  calculated  from  the  grid  array 
data.  The  parameter  results  were  calculated  for  a specific  water  ele- 
vation by  calculating  the  area  of  the  bottom  surface  within  each  range 
class  for  that  water  elevation.  The  computer  calculation  procedure 
used  for  each  water  elevation  was  follows.  A series  of  six  horizontal 
planes  were  constructed  at  2,  4,  6,  10,  15,  and  20  ft  below  the  water 
surface . Each  grid  square  within  the  side  channel  was  broken  into  two 
triangles  by  placing  a diagonal  on  the  square,  thus  forming  two  trian- 
gular planes . Each  triangle  was  then  further  broken  into  10  parallelo- 
grams and  each  parallelogram  was  assigned  an  average  elevation.  The 
average  elevation  of  each  parallelogram  was  then  examined  to  determine 
whether  it  was  above  water  level,  and  therefore  ignored,  or  below  water 
level.  If  it  was  below  water  level,  it  was  examined  to  determine  which 
depth-class  range  it  belonged  to,  and  the  area  of  that  parallelogram  was 
added  to  the  accumulated  total  area  for  that  class  range.  The  operation 
is  expressed  mathematically  by  the  following  equation: 


1* 

-12 


d2A.  . 
ij 


where 


B.  = the  bottom  area  in  range  class  j when  the  water  elevation 
^ is  at  the  k^h  level 

N = the  number  of  triangles  in  the  grid  array 
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class  range 


d = the  grid  spacing 

A.  = the  fraction  of  the  i ° triangle  in  the 

75.  The  volume  of  information  produced  by  performing  this  calcu- 
lation is  large,  and  the  parameter  value  trends  are  not  immediately 
apparent  when  viewing  a table  of  values.  To  partially  remedy  that 
situation,  the  data  were  also  expressed  in  a normalized  form  in  which  the 
data  for  each  water  elevation  were  normalized  to  the  total  bottom  surface 
area  underwater.  The  mathematical  procudure  followed  is  shown  below, 
and  both  the  absolute  and  normalized  forms  of  the  data  are  presented  in 
the  results. 


, 100. OB. 

J T 


B 


i= i 


where 


= the  percentage  of  the  total  bottom  area  underwater  in  the 
class  range  j when  the  water  is  at  the  level 


Parameter  k:  Water  cross-sectional 

area  as  a function  of  water  eleva- 
tion  at  selected  sampling  locations 

76.  The  water  cross-sectional  area  as  a function  of  water  eleva- 
tion was  calculated  from  the  elevation  grid  array  data.  Calculations 
were  performed  only  at  specific  locations  (stations).  The  river  miles 
of  the  locations  are  listed  in  table  1 and  the  stations  are  shown  in  the 
contour  maps  in  Volume  II. 

77*  The  calculational  procedure  was  the  same  for  all  stations.  A 
vertical  plane  was  placed  perpendicular  to  the  side-channel  center  line 
at  the  station,  and  the  profile  formed  by  the  intersection  of  that  plane 
and  the  side-channel  bottom  was  calculated.  Horizontal  planes  were  then 
placed  at  2-ft  intervals  through  that  profile  and  the  area  bounded  by 
each  plane  and  the  bottom  of  the  profile  was  calculated.  The  trapezoidal 
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rule  was  used  for  the  area  calculation.  The  mathematical  equation  is  as 
follows : 


the  cross-sectional  area  of  the  water  in  the  side  channel  when 
the  water  is  at  the  k^h  level 

the  first  point  in  the  string  of  coordinates  defining  the  pro- 
file that  is  above  water  as  the  profile  advances  upward  toward 
the  high  bank  from  the  side-channel  bottom 

the  last  point  in  the  string  of  coordinates  defining  the  pro- 
file that  is  above  water  as  the  profile  advances  down  toward 
the  water  from  the  high  bank 

the  coordinates  of  profile  in  a plane  that  cuts  the  side 
channel  perpendicular  to  the  center  line 


Closure  of  the  curve  during  the  areal  calculation  was  satisfied  by  the 
condition  ^4+1  ,ZM+1  ^ = The  calcula'ti°n  °f  the  water  cross 

section  over  only  the  water-covered  region  of  the  profile  was  satisfied 
by  applying  the  conditions  Z.  <_  E or_  Z < i during  calculation. 


78.  The  results  of  all  calculations  are  presented  in  tables  2-19 
and  plates  1-18.  Each  table  contains  the  total  results  for  a side 
channel.  Each  plate  is  a graphic  presentation  of  a portion  of  the 
tabular  data  for  a side  channel.  The  tabular  data  not  represented  in 
plates  l-l8  are  the  length  along  thalweg,  the  average  width  between  high 
banks,  and  the  side-channel  bottom-surface  area  underwater  as  a function 
of  water  elevation  and  water  depth.  The  former  two  parameters  are  one- 
dimensional and  cannot  be  represented  graphically.  The  lattermost 
parameter  is  three-dimensional  and  could  have  been  represented 


graphically,  but  was  not,  since  no  advantage  could  be  gained  because  of 
the  difficulty  in  interpreting  such  a graphic. 

79.  The  tabular  and  graphic  results  are  arranged  by  side  channels 
as  shown  below: 


Table 


Plate 


Side  Channel 

No. 

No. 

Jefferson  Barracks 

2 

1 

Calico 

3 

2 

Osborne 

hA 

3A 

Osborne,  Pool  1 

kB 

3B 

Osborne , Pool  2 

hC 

3C 

Osborne,  Pool  3 

LD 

3D 

Osborne,  Pool  U 

Ue 

3D 

Osborne , Pool  5 

Uf 

3E 

Harlow 

5 

k 

Fort  Chartres 

6a 

5A 

Fort  Chartres , Pool  1 

6b 

5B 

Fort  Chartres,  Pool  2 

6c 

5B 

Fort  Chartres,  Pool  3 

6d 

5C 

Moro 

7 

6 

Kaskaskia 

8 ' 

7 

Crains 

9 

8 

Liberty 

10A 

9A 

Liberty,  Pool  1 

10B 

9B 

Liberty,  Pool  2 

IOC 

9B 

Liberty,  Pool  3 

10D 

9C 

Liberty,  Pool  h 

10E 

9C 

Liberty,  Pool  5 

10F 

9D 

Jones 

11 

10 

Picayune 

12 

11 

Cape  Bend 

13 

12 

Sante  Fe 

lU 

13 

Billings 

15 

lh 

Buffalo 

16 

15 

Browns 

17 

16 

Thompson 

18 

17 

Sister 

19 

18 

80.  Values  for  all  parameters  are  given  in  the  tables  for  both  the 


total  side  channel  and  also  for  the  pools  (except  parameter  k)  of  three 
selected  (Osborne,  Fort  Chartres,  and  Liberty)  side  channels.  The 
parameters  that  are  functions  of  water  elevation  are  tabulated  both 
according  to  mean  sea  level  and  according  to  St.  Louis  gage  equivalent 
water  elevations  at  the  side  channel.  Each  table  is  formatted  in  an 
identical  manner.  Results  for  parameters  a and  b_  (center-line  length 
and  average  width  between  high  banks)  appear  at  the  top  of  the  table, 
those  for  parameters  c_  through  i_  (volume,  shoreline  length,  water  surface 
area,  shoreline  development,  rate  of  change  of  water  surface  area  with 
respect  to  water  elevation,  ratio  of  water  surface  area  to  volume,  and 
ratio  of  shoreline  length  to  water  surface  area)  appear  in  the  upper 
box,  those  for  parameter  (bottom  surface  area  underwater  as  a function 
of  water  elevation  and  water  depth)  in  the  lower  box,  and  those  for 
parameter  k_  (water  cross-sectional  area  as  a function  of  water  elevation 
at  selected  sampling  locations)  in  the  boxes  to  the  right  in  each  table. 

8l.  The  graphic  representations  of  the  data,  as  seen  in 
plates  1-18,  were  produced  to  provide  data  for  those  parameters  in  a form 
immediately  useful  in  the  first  step  required  to  correlate  geometric 
characteristics  with  wildlife  populations,  and  to  present  them  in  a 
particularly  simple  and  readily  interpretable  form.  The  intention  was 
to  provide  the  data  so  that  trends  in  parameters  as  a function  of  water 
elevation  could  be  easily  seen,  and  the  trends  for  one  side  channel 
could  be  compared  with  those  of  any  other.  Since  each  side  channel  had 
a different  size  and  shape,  the  absolute  values  of  the  parameters  were 
different  for  each.  To  provide  a uniform  means  of  displaying  the 
selected  data  that  is  consistent  with  the  purpose,  all  selected  data  are 
displayed  with  a normalized  format.  Normalization  of  data  for  each  para- 
meter of  each  side  channel  was  performed  relative  to  the  maximum  value  for 
that  parameter  for  that  side  channel.  The  mathematical  procedure  was  as 
follows : 


= 100 


\(MAX) 


where 


of  maximum. 


PART  IV:  CONCLUSIONS  AND  RECOMMENDATIONS 


\ 


Conclusions 


82.  A procedure  for  quantifying  geometric  characteristics  of 
water  bodies  was  successfully  developed  and  used  in  the  study  of  18 
Middle  Mississippi  River  side  channels.  Values  for  those  parameters 
selected  to  describe  the  geometric  characteristics  quantitatively  were 
calculated  using  assorted  available  data.  Since  the  statistical  sample 
size  used  to  calculate  the  value  of  any  of  the  parameters  ranged  from 
3,000  to  10,000,  the  calculated  values  of  the  parameters  can  be  con- 
sidered more  accurate  than  the  values  of  parameters  calculated  for 
chemical  or  biological  characteristics  (the  sample  sizes  available  for 
those  characteristics  are  simply  much  smaller) . The  primary  purpose  of 
this  study  was  to  provide  data  for  comparing  the  geometric,  chemical, 
and  biological  elements  of  a side  channel.  The  accuracy  of  the  data 
calculated  in  this  study  is  sufficient  for  that  purpose. 

Recommendations 


83.  No  further  calculations,  either  of  the  same  parameters  in 
additional  side  channels  or  more  parameters  for  the  18  side  channels, 
should  be  undertaken  until  the  results  of  this  study  are  used  for  the 
intended  purpose.  The  correlation  of  the  calculated  parameter  data  with 
wildlife  populations  should  be  performed  and  studied  to  determine  which 
of  the  parameters  are  best  indicators  of,  for  example,  total  fish  popula- 
tions and  densities.  The  first  and  most  basic  step  in  that  correlation 
should  be  the  grouping  of  side  channels  according  to  similarity  of  param- 
eter function  shapes  using  plates  1-18.  Since  more  advanced  classical 
correlation  coefficient  calculations  are  neither  straightforward  nor 
easily  understood  for  two-  and  three-dimensional  parameters,  attempts 
should  be  made  to  show  correlations  on  gross  features  (e.g.  total  fish 
population  density  irrespective  of  species  against  shoreline  development) 
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to  aid  in  understanding  the  data  and  possibly  to  uncover  simple,  strong 
relations,  if  there  are  any. 

8U.  The  data  can  also  be  rearranged  to  yield  nev  insights.  Hie 
monthly  river-stage  records  in  the  immediate  locale  of  each  side  channel 
can  be  used  in  conjunction  with  the  elevations  of  controlling  structures 
to  produce  estimates  for  the  parameters  as  a function  of  month  rather 
than  water  level.  Parameter  values  at  approximate  fish  spawning, 
juvenile,  and  mature  stages,  and  duck  migration  times  could  be  used  to 
classify  the  side  channels  according  to  similarity. 

85.  Portions  of  the  data  could  be  used  for  other  purposes.  New 
fathometer  and  aerial  photo  data  could  be  used  to  perform  calculations 
for  selected  side  channels  to  study  the  change  in  shape,  particularly 
changes  possibly  due  to  the  1972-73  flood.  Contour  maps  for  different 
times  could  synthesize  the  changes  for  ease  of  interpretation.  Even 
significant  changes  in  a side-channel  contour  map  would  not  necessarily 
imply  significant  changes  in  the  parameter  values. 
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Table  1 


Side  Channels  Included  in  Study 


Name 

No. 

River  Mile 
Extreme s* 

River 

1 

Mile  Station 
2 

Locations 

3 

Jefferson  Barracks 

1 

168. 8-166. 5L 

166.6 

167.5 

168.6 

Calico 

2 

148. 3-147. 4L 

147.5 

147.9 

148.1 

Osborne 

3 

l46.3-l44.3L 

144.6 

145-3 

145.8 

Harlow 

4 

l43.0-l4l.7R 

142.0 

142.3 

142.7 

Salt  Lake 

5 

141. 5-136. 8L 

137.8 

138.8 

139.9 

Fort  Chartres 

6 

134. 3-132. 3L 

132.8 

133-1 

134.1 

Establishment 

7 

132. 5-130 .1R 

130.9 

131.5 

132.3 

Moro 

8 

122. 6-120. 1L 

120.0 

121.7 

122.3 

Kaskaskia 

9 

118. 0-115- 8R 

116.5 

117.2 

117.4 

Crains 

10 

105. 6-104. 4r 

104.5 

105.1 

105.4 

Liberty 

11 

102. 8-99. 9L 

100.4 

101.3 

102.0 

Jones 

12 

98. 3-95. OR 

95.6 

96.4 

98.2 

Grand  Tower 

13 

78. 5-77. 8R 

-- 

-- 

-- 

Crawford 

14 

73. 7-71. 6l 

72.2 

72.6 

73-6 

Schenimann 

15 

62. 5-57. 1R 

57.7 

59-1 

61.7 

Picayune 

16 

60. 6-54. 8L 

54.9 

58.0 

60.1 

Cape  Bend 

17 

51. 4-47. 6l 

48.0 

48.8 

50-7 

Santa  Fe 

18 

4o.4-35.4l 

35.6 

-- 

37-3 

Billings 

19 

34. O-32. 7R 

33-0 

33.5 

33-9 

Buffalo 

20 

26.8-"4.7R 

25.O 

25.5 

26.0 

Browns 

21 

24. 7-21. 7L 

22.1 

22.8 

24.2 

Thompson 

22 

18. 7-15. 3R 

15.9 

17.7 

18.6 

Sister 

23 

l4.4-U.9R 

12.5 

13.0 

14.0 

Boston 

24 

10.4-7.7L 

8.3 

8.8 

10.1 

Angelo 

25 

5-0-1. 2L 

2.6 

4.0 

4.8 

* R and  L indicate 

that 

the  side  channel 

is  on  the  right  or  ' 

left  of 

the  river,  respectively,  when  facing  downstream. 


Results  of  Calculations  for  Side  Channel  Osborne 


Table  4 (Continued) 


Table 


Table  4 (Continued) 

E.  Side  Channel  Osborne,  Pool  4 

s I OE  CHANNEL  LENGTH  * 0.38  HUES  * 0.61  K ILOHETENS 


Table  4 (Concluded.) 


Result,  of  Calculations  for  Side  Channel  Fort  Chartres 


Table  6 (Continued) 

Side  Channel  Fort  Chartres,  Pool  1 


Table  6 (Continued) 

Side  Channel  Fort  Chartres,  Pool  2 


calculated  depth  class  ranges 


CALCULATED  DEPTH  CLASS  RANGES 


Results  01'  Calculations  for 


CALCULATED  DEPTH  CLASS  RANGES  CALCULATIONS  OF  PROFILE  CROSS  SECTION 


CALCULATED  DEPTH  CLASS  RANGES  CALCULATIONS  OF  PROFILE  CROSS  SECTION 


Results  of  Calculations  for  Side  Channel  Liberty 


Side  Channel:  Liberty,  Pool  2 


Table  10  (Continued) 


Table  10  (Continued) 

Side  Channel:  Liberty,  Pool  4 


Table  10  (Concluded) 


Results  of  Calculations  for 


CALCULATED  DEPTH  CLASS  RANGES  CALCULATIONS  OF  PROFILE  CROSS  SECTION 


Results  of  Calculations  for 


CALCULATED  DEPTH  CLASS  RANGES  CALCULATIONS  OF  PROFILE  CROSS  SECTION 


Results  of  Calculations  for 


Results  of  Calculations  for 


CALCULATED  DEPTH  CLASS  RANGES  CALCULATIONS  OF  PROFILE  CROSS  SECTION 


Results  of  Calculations  for 


CALCULATED  DEPTH  CLASS  RANGES  CALCULATIONS  OF  PROFILE  CROSS  SECTION 


Side  Channel  Browns 


CALCULI  TIP  DEPTH  CLASS  ffAWGES  CALCULATIONS  OF  PROFILE  CROSS  SECTION 


Results  of  Calculations  for 


25  13  5 235  36  0 

RIVER  STAGE,  FT 


PARAMETER  VALUES  FOR 
ST.  LOUIS  GAGE 
EQUIVALENT  RIVER  STAGE 
FOR  MISSISSIPPI  RIVER 

SIDE  CHANNEL 
JEFFERSON  BARRACKS 


2 5 13  5 23  5 360 

RIVER  STAGE.  FT 


PLATE  1 


RIVER  STAGE,  FT 


PARAMETER  VALUES  FOR 
ST.  LOUIS  GAGE 
EQUIVALENT  RIVER  STAGE 
FOR  MISSISSIPPI  RIVER 
SIDE  CHANNEL  CALICO 


30  14  0 25  0 

RIVER  STAGE  . FT 


PLATE  2 


<-3.5  <-3  5 <-35  -15  10.0  21.0 

RIVER  STAGE,  FT 


PARAMETER  VALUES  FOR 
ST  LOUIS  GAGE 
EQUIVALENT  RIVER  STAGE 
FOR  MISSISSIPPI  RIVER 

SIDE  CHANNEL  OSBORNE 


0 -i r 1 r 1 1 

<-35  <-3  5 <-3  5 -15  10.0  210 

RIVER  STAGE.  FT 


PLA I h 3A 


PARAMETER  VALUES  FOR 
ST  LOUIS  GAGE 
EQUIVALENT  RIVER  STAGE 
FOR  MISSISSIPPI  RIVER 

SIDE  CHANNEL  OSBORNE 
POOL  I 


100  210 
RIVER  STAGE . ET 


PLATE  3B 


<-35  <-35  <-3  5 -15  100 

RIVER  STAGE.  FT 


PARAMETER  VALUES  FOR 
ST  LOUIS  GAGE 
EQUIVALENT  RIVER  STAGE 
FOR  MISSISSIPPI  RIVER 

SIDE  CHANNEL  OSBORNE 
POOL  2 


PLATE  3C 
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